Journal of Biomechanical Engineering

Repetitive low-intensity blast (LIB) exposure has been identified as a probable cause of mild blast-induced traumatic brain injury (mbTBI) and a chronic injury risk to U.S. military personnel. However, the human brains biomechanical response to this loading regime remains poorly characterized. Using blast-exposure-validated 3D models of human anatomy, we simulated the intracranial tissue response to blast pressure typically experienced by Warfighters during weapons training. Two scenarios were evaluated, a single-dose exposure and a repetitive-dose exposure, to study intracranial pressure (ICP), shear strains, and spectral content. Ansys LS-DYNA was used to generate planar blast waves with peak overpressures of 4-90 kPa and positive phase durations of 2.2-10 ms. Single exposures produced ICP ranging from 4.7-112.7 kPa, dependent on dose and positive phase duration. Under repetitive LIB exposure, peak ICP increased by 8-26% relative to single exposures, with an increase of high-frequency components (>2 kHz). These results demonstrate that LIB can produce measurable intracranial responses that are amplified through repetition, producing pronounced spectral content and elevated pressures despite low strain levels. This study underscores the need to further investigate cumulative dose effects and the value of computational approaches to clarify hypothesized mbTBI mechanisms in operationally relevant conditions.

Adolescent idiopathic scoliosis (AIS) alters spinal geometry and may influence the biomechanical response of the spine during functional postures. However, posture-dependent changes in spinal loading and paraspinal muscle forces in AIS remain poorly understood. This study investigated the effects of trunk posture on intervertebral loading and paraspinal muscle forces using a subject-specific musculoskeletal model of an adolescent with AIS. The spinal deformity was reconstructed from biplanar radiographs and incorporated into a full-body musculoskeletal model. Flexion, extension, lateral bending, and axial rotation were simulated at three incremental magnitudes, with motion distributed across the thoracolumbar spine. Intervertebral compressive and lateral forces around the curve apex and forces in the erector spinae (ES) and multifidus (MF) muscles were evaluated. Trunk flexion produced the greatest compressive loading, reaching 337 N at the curve apex and 372 N two levels below the apex at 30{degrees} flexion. Lateral bending produced pronounced direction-dependent loading: concave-side bending increased lateral forces, whereas convex-side bending increased compressive forces. Axial rotation produced similar but smaller direction-dependent changes. Paraspinal muscle forces were consistently asymmetric, with concave-side dominance of the ES and convex-side dominance of the MF. Flexion and convex-sided movements generally produced greater muscle imbalance, while increasing posture magnitude amplified spinal loading and muscle forces. These findings demonstrate that trunk posture, movement direction, and magnitude substantially influence the biomechanical environment of the scoliotic spine and should be considered when evaluating spinal mechanics in AIS.

Ultrasound is commonly used to assess structural changes in symptomatic Achilles tendons, but quantitative biomechanical metrics for progressive tendon deterioration remain limited. The goal of this study was to develop and validate an automated ultrasound tracking algorithm for regional tendon deformation and evaluate strain progression in survived and ruptured tendons during fatigue loading. We hypothesized that maximum strain, average strain, and strain heterogeneity would exhibit different trajectories between groups. Ten cadaveric Achilles tendons underwent cyclic loading with stress tests every 500 cycles until rupture or 150,000 cycles. Ultrasound images acquired during stress tests were analyzed using an automated tracking algorithm to generate spatially resolved regional strain fields. Ultrasound-derived bulk strain was highly correlated with actuator-derived strain in survived (R^2 = 0.968 +/- 0.017) and ruptured tendons (R^2 = 0.972 +/- 0.014). Maximum and average longitudinal strains progressively diverged between groups across fatigue life (Group x FatigueLife: p = 0.003 and p < 0.0001, respectively). During the first 10,000 cycles, average strain decreased in survived tendons ({beta} = -0.0268%, p = 0.0215) but not ruptured tendons ({beta} = 0.0147%, p = 0.1197), with a significant Group x Cycle interaction (p = 0.0061). This study demonstrates that the algorithm quantified Achilles tendon deformation with high fidelity and enabled spatially resolved strain assessment throughout fatigue loading. Maximum and average strain followed different trajectories between groups, whereas strain heterogeneity did not. Early differences in tendon biomechanics suggest that regional strain behavior may change before pronounced differences in absolute magnitude develop.

This study presents a physiologically detailed biomechanical model of the mouse distal forelimb that incorporates intrinsic musculature, tendon routing, and digit-level skeletal anatomy, features simplified or omitted in existing musculoskeletal models. Using high-resolution anatomical reconstruction and computational modeling, we created a physiological representation of the wrist and digits capable of simulating complex forelimb movements. The model enables simulation of coordinated distal forelimb movement and digit-level muscle behavior during grasping-related tasks. Simulations were performed for multiple tasks, including grasping, grasping with supination, wrist flexion, and digit I flexion, with analysis focused on the grasping task due to its integration of both intrinsic and extrinsic musculature. Model performance was evaluated through comparisons of marker trajectories between torque-driven reference motion and muscle-driven simulations, temporal shuffle control, and comparisons between experimentally recorded electromyography (EMG) activity and model-predicted muscle excitation profiles. The model successfully reproduced coordinated distal forelimb kinematics, demonstrated strong agreement between torque-driven and muscle-driven simulation approaches, and generated physiologically plausible muscle excitation patterns consistent with experimentally observed EMG activity during grasping-related movement. These findings establish the model as a framework for studying fine motor control, neuromuscular coordination, and movement-related impairments in mice while providing a foundation for future investigation of neurological disorders and their underlying biomechanical mechanisms.

Successful instrumented fusion of the lumbar spine is a complex surgical challenge, with positive patient outcomes dependent on careful surgical planning. Material selection is of critical importance to a mechanical construct supporting successful spinal fusion. Therefore, the aims of this study were to (a) evaluate the potential clinical use of finite element analysis (FEA) and (b) conduct a retrospective mechanical analysis of different implant materials in patients having undergone spinal fusion using FEA. Our methodology involved segmenting the spine from post-operative computed tomography (CT) image data from patients with previous spinal fusion. FEA models representing post-surgery cases were developed and different biomechanical loading conditions such as compression, flexion, bending and extension whilst testing pedicle screws of different materials were simulated. Patient specific finite element models were created, and biomechanical analysis were completed for all three patients. Polyetheretherketone (PEEK) constructs typically demonstrated lower peak implant stress when compared to titanium constructs for all spinal fusion levels. Furthermore, increasing the spinal fusion level resulted in significant differences in the maximum von Mises stress within both the bone and the instrumentation, whereas the 2-level fusion exhibited comparable stress levels in the bone irrespective of the instrumentation material. This pilot explores the potential of FEA as a clinical tool for assessing device and bone stresses. In our cohort, different materials can influence the stresses in both the instrumentation and the instrumented vertebrae, suggesting FEA can be useful pre- operative tool with regards to instrument selection and post-operatively to assess instrumentation and bone stresses. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=74 SRC="FIGDIR/small/742715v1_ufig1.gif" ALT="Figure 1"> View larger version (33K): org.highwire.dtl.DTLVardef@55959aorg.highwire.dtl.DTLVardef@d0b9d6org.highwire.dtl.DTLVardef@158c348org.highwire.dtl.DTLVardef@7ce828_HPS_FORMAT_FIGEXP M_FIG C_FIG

Predicting wear in artificial joints requires integrating joint dynamics, contact mechanics and progressive surface evolution, yet these processes are often treated separately. In total hip arthroplasty (THA), finite-element approaches remain the reference standard, but they are computationally demanding and usually rely on boundary conditions from independent musculoskeletal (MSK) models, hindering consistent coupling and feedback between wear progression and movement dynamics. As a single-subject proof of concept, we present a computational framework that embeds wear estimation within forward MSK simulations through OpenSim-MATLAB integration. Contact variables are computed using an elastic-foundation formulation, and wear is updated through the Archard law, enabling cyclic prediction of contact mechanics and surface evolution within a single workflow at practical computational cost. The framework was evaluated in one subject with right THA during five activities of daily living and numerically benchmarked against finite-element simulations. A long-term walking analysis of 4 million cycles was also performed to assess geometry updating. Across tasks, peak contact pressures remained within 7% of finite-element predictions. Linear wear depth and volumetric loss showed maximum deviations of 16% and 13%, respectively. Accounting for progressive geometry changes yielded a maximum wear depth about 31% lower than linear extrapolation. These preliminary results support the framework's computational feasibility and numerical consistency for the tested case; nevertheless, multi-subject evaluation is required before broader predictive or clinical use.

Background: Ascending thoracic aortic aneurysm (ATAA) is characterized by progressive biomechanical remodeling that is incompletely captured by maximum aortic diameter. The relationships between age, sex, aortic valve phenotype (BAV vs TAV), and aneurysmal mechanics remains incompletely understood. Methods: Ascending aortic specimens from 54 patients undergoing elective ATAA repair were evaluated. Biaxial testing quantified Green strain and local stiffness at membrane tensions of 60 and 120 N/m. Opening angle was used to assess circumferential residual stress and stretch. Associations with age were evaluated using Pearson correlations, and differences between sex and valve phenotype were also explored. Multivariable linear regression models including age, sex, and valve phenotype were used to evaluate their adjusted associations with each mechanical outcome. Results: Increasing age was associated with increased local stiffness and reduced Green strain and remained the strongest determinant of biaxial mechanical properties after adjustment. Compared with tricuspid aortic valve specimens, BAV specimens exhibited lower Green strain together with lower local stiffness. Females demonstrated greater Green strain in unadjusted analyses; however, after adjustment, sex remained associated only with longitudinal Green strain at 60 N/m. None of the evaluated clinical characteristics consistently explained opening angle or residual stretch. Conclusions: Age was the predominant clinical determinant of biaxial mechanical remodeling in ATAA, whereas valve phenotype demonstrated selective associations and sex contributed comparatively little after adjustment. Circumferential residual stress appeared largely independent of these clinical characteristics, suggesting that residual stress and biaxial mechanics capture complementary aspects of aneurysmal remodeling.

Introduction: The femoropopliteal artery undergoes substantial axial compression, bending, torsion, and pinching during limb flexion. When an implanted stent does not accommodate these deformations, repetitive biomechanical mismatch may alter the vascular response to implantation and contribute to restenosis. We tested whether designing a peripheral stent to better accommodate physiologic limb-flexion deformation reduces the subsequent inflammatory and neointimal response. Methods: A nitinol stent was developed using multi-objective computational optimization to maintain arterial apposition while reducing wall stress and luminal pinching during physiologic deformation. Sixteen Yucatan minipigs were studied in independent 3- and 6-month cohorts (n=8 each). Each animal received the optimized stent in one limb and a commercial stent with representative biomechanical behavior contralaterally. CTA quantified in vivo limb-flexion deformation, while CT and histology assessed stenosis, neointimal accumulation, inflammatory-cell burden, SMA-positive tissue, and extracellular-matrix composition. Results: The optimized stent accommodated greater axial shortening than the commercial comparator at both 3 months (11.2{+/-}2.7% vs. 6.6{+/-}2.1%, p=0.003) and 6 months (13.7{+/-}4.4% vs. 6.8{+/-}2.3%, p=0.004). Intimal hyperplasia volume was lower with the optimized stent at 3 months (227.5{+/-}85.2 vs. 350.0{+/-}130.5 mm3, p=0.002) and 6 months (183.6{+/-}46.2 vs. 261.2{+/-}60.6 mm3, p=0.001), with significantly lower maximum, mean, and median stenosis at both timepoints. At 3 months, macrophage and total inflammatory-cell counts were 51% and 45% lower, respectively, with the optimized stent. Histologic tissue area was 28% lower at 3 months and 12% lower at 6 months. SMA-positive tissue fraction was higher with the optimized stent at both timepoints, whereas elastin, glycosaminoglycan, and collagen fractions were similar between devices. Discussion: Improved accommodation of limb-flexion deformation was associated with a smaller inflammatory and neointimal response to stent implantation through 6 months, while overall neointimal matrix composition remained largely unchanged. These findings provide in vivo evidence that biomechanical compatibility with the dynamic lower-limb environment can influence vascular healing and should be considered a mechanistic design criterion for peripheral artery stents.

Purpose: Fibrotic remodeling of the uterus, associated with aging, disease, and environmental exposures, alters collagen organization and tissue stiffness, yet how these changes influence organ-level mechanical behavior remains poorly understood. Glutaraldehyde (GA)-induced collagen crosslinking was used as a controlled surrogate for fibrotic remodeling to determine whether image-informed inverse finite element analysis (iFEA), combined with inflation testing and micro-computed tomography (microCT), could detect and quantify the resulting changes in uterine constitutive behavior. Methods: Murine uteri (n = 6 untreated, n = 6 GA-crosslinked) underwent volume-controlled balloon inflation with simultaneous microCT imaging to quantify deformation of the inner and outer wall boundaries for iFEA. Specimen-specific Gasser-Ogden-Holzapfel (GOH) finite element models were optimized by adjusting model parameters to reproduce experimentally measured wall contours throughout inflation. Model performance was evaluated using contour root mean square error (RMSE), and parameter identifiability was assessed through sensitivity analyses. Results: GA treatment significantly increased inflation work, linear stiffness, and maximum inflation resistance (p < 0.001). The iFEA framework accurately reproduced experimental deformation (RMSE < 3%) and revealed significant increases in the estimated GOH parameters C10 (9.2-fold), k1 (2.0-fold), and k2 (2.7-fold), consistent with increased effective tissue stiffness and a shift toward earlier collagen fiber recruitment. Sensitivity analyses demonstrated unique, well-defined minima for all parameter combinations. Conclusion: Image-informed iFEA provides a quantitative framework for relating collagen remodeling to organ-level uterine mechanics through specimen-specific constitutive parameter estimation. This approach establishes a foundation for investigating the mechanical consequences of uterine fibrosis and other remodeling processes.

Healing of the gastrointestinal (GI) tract following surgical interventions is governed by coupled mechanobiological processes and may lead to severe complications such as perforation or stricture formation. A first comprehensive mechanobiological modeling framework of GI healing is herein proposed. The framework addresses two clinically relevant applications of bioprinting: endoscopic treatment of GI tissue defects and bioprinting-assisted surgical anastomosis. Tissue remodeling is based upon a hyperelastic material model with fiber dispersion coupling the spatiotemporal reaction diffusion dynamics of fibroblasts and cytokines with collagen deposition, reorientation, and plastic deformation. An active stress formulation drives fibroblasts and myofibroblast contraction while a multiplicative decomposition approach is adopted for collagen remodeling as a permanent change. The governing equations are discretised in time with an implicit Euler scheme, finite element in space and solved by a staggered algorithm implemented in the open-source platform FEniCS. Numerical simulations show that the initial collagen fibers content and orientation strongly influence fibroblast infiltration and scar contracture formation. The study quantifies the advantage of having a collagen bioprinted architecture aligned with tissue fibers, to reduce structural remodeling distortion. The model also predicts a strong correlation between collagen deposition and luminal burst pressure during the recovery of mechanical integrity in GI anastomosis.

Finite element (FE) head models are virtual tools to study brain biomechanics and their predictions must be numerically convergent. Previous convergence studies focused on the influence of mesh size, but the potential effect of element formulation on model convergence was often ignored. To address this, one original model with brain mesh size as 6.4 {+/-} 1.9 mm was modified to generate three derivatives with the same mesh topology but different element sizes, i.e., a coarse model (mesh size: 12.2 {+/-} 3.9 mm), a medium model (mesh size: 3.2 {+/-} 1.0 mm), and a fine model (mesh size: 1.6 {+/-} 0.5 mm). Three commonly used element formulations, i.e., reduced integration, selectively reduced (S/R) integration, and full integration, were implemented to the brain elements. These models were subjected to rotational loadings along the axial, coronal, and sagittal axes, respectively. The maximum relative displacement at representative sites and 95th percentile maximum principal strain at the whole brain level were used to evaluate mesh convergency. The results showed that the S/R integration yielded a 5% difference between the original and medium meshes, while the reduced and full integration revealed a difference over 5% even between the medium and fine meshes. This study verified that the mesh convergence of FE brain models is affected by the choice of element formulation and the S/R integration contributes to the fastest convergence behavior than the reduced and full integrations. It provided practical information on how to develop numerically convergent and computationally efficient FE brain models. HighlightsO_LIThis study verifies that the choice of element formulation affects the mesh convergence behavior of finite element brain models C_LIO_LIThis study finds the selectively reduced integration yields the fastest convergence behavior than the reduced and full integration C_LIO_LIThis study provides practical guidance on the choice of mesh density and element formulation on how to develop numerically convergent and computationally efficient finite element brain models. C_LI

The primary source of strength in ligaments and tendons are hierarchically organized collagen fibers. These fibers largely do not regenerate after injury, with repair, nor in engineered replacements, limiting treatment options. Previously, we developed a culture system which guides ACL fibroblasts in high-density collagen gels to form native-size hierarchical fibers over 6 weeks, and demonstrated that intermittent cyclic stretch further improves maturation. However, additional maturation is needed for clinical relevance. Interestingly, we found cyclic load affected cells differentially depending on the degree of organization, with 10% cyclic strain driving early improvements in unorganized gels and 5% strain being more beneficial later in culture once cells were on aligned fibers. Here, we explored whether a stepped cyclic load, that increased or decreased in strain magnitude as collagen fibers developed, further improved maturation. We hypothesized that progressively decreasing cyclic strain as organization increases would drive cells to produce more mature hierarchical fibers, resulting in stronger replacements. Controls had intermittent cyclic stretch at 0, 5, 7, or 10% strain throughout culture, while stepped load constructs were cyclically loaded with a strain that increased or decreased by 2-3% every 2 weeks as constructs matured. Contrary to our hypothesis, neither decreasing nor increasing load led to further tissue maturation. We hypothesize stepped cyclic load may disrupt cellular tensional homeostasis, leading to repeated remodeling of collagen and shifted proteoglycan accumulation. This study provides insight into how stepped cyclic loading affects hierarchical fiber formation and maturation, which will help to engineer stronger replacements and better rehabilitation protocols.

ObjectiveThis study investigates knee joint biomechanics during lunges under varying tibial angles and external loads using musculoskeletal modeling and finite element analysis. The goal is to provide a biomechanical basis for understanding knee loading patterns and optimizing sports training and rehabilitation. MethodsTwenty-six healthy young men performed lunges under tibial inclination angles relative to the ground (60{degrees} and 90{degrees}) and two external load conditions (bodyweight and an additional 98 N external load). Kinematic and kinetic data were captured using motion capture and force plates. Musculoskeletal models were used to estimate joint moments, range of motion, and stiffness, with data analyzed using two-way repeated-measures ANOVA. Finite element analysis was performed at 90{degrees} tibial angle to evaluate tissue stress and displacement. ResultsThe joint moment at a 60{degrees} tibial angle was much higher than at a 90{degrees}. External load showed significant effects on knee stiffness, with lower rotational stiffness in the horizontal plane (P < 0.001) and lower coronal plane stiffness at 90{degrees} (P = 0.012) under the 98 N external-load condition, indicating reduced resistance to angular displacement in these planes. Under the 90{degrees} tibial-angle condition with external load, peak stress and displacement were concentrated in the posterior horn of the meniscus, with a maximum displacement of 3.12 mm. ConclusionThe anterior tilt of the tibia increased sagittal-plane knee loading, while external load mainly reduced joint stiffness in the coronal and horizontal planes. Under the 90{degrees} loaded condition, the concentration of stress and displacement in the posterior horn of the meniscus suggests a mechanically unfavorable loading pattern rather than direct evidence of injury risk. These findings may provide useful biomechanical information for load management during lunge-based training and rehabilitation.

To correct lower limb length discrepancies, orthopedic surgeons often perform limb lengthening procedures. These procedures involve directly lengthening a primary bone (e.g., femur), which also stretches the muscles and tendons spanning the elongation site. Limb lengthening procedures improve lower-limb length symmetry, often at the expense of chronic joint-level decrements. For example, knee joint range of motion and strength does not always return to pre-surgery values following femur lengthening. These joint-level decrements governed by altered muscle force-length characteristics, but it is unestablished how human muscles remodel during surgical limb lengthening. Here, we characterized bilateral human vastus lateralis fascicle force-length profiles following femur lengthening. At 2.2 {+/-} 1.7 years (Avg {+/-} SE) post-initial surgery, the lengthened limb's vastus lateralis fascicles were longer, effectively stiffer, and weaker than the contralateral limb. The lengthened limb's vastus lateralis fascicles were 38 mm longer (Avg) than that of the contralateral limb, even though the lengthened limb's femur was only 31 mm longer (Avg). The onset of vastus lateralis fascicle passive force production began at 22% shorter (Avg) operating lengths in lengthened versus contralateral limb muscles. Additionally, vastus lateralis peak active force production was 34% less (Avg) in the lengthened versus contralateral limbs, and in-series quadriceps tendon stiffness was not stiffer in lengthened limbs. These data suggest that active muscle tissues lengthen more than passive muscle tissues following limb lengthening, effectively creating stiffer muscles. Years after limb lengthening, patient muscles exhibit numerous differences in their force-length properties versus that of natural growth, likely contributing to impaired joint-function.

Despite lacking a dedicated structural function, brain tissue is exposed to severe mechanical loading during accidents and trauma, with potentially permanent or lethal consequences. Therefore, characterizing its mechanical behavior assumes particular importance for designing protective equipment, planning surgeries, and training neurosurgeons, tasks that often rely on physical brain replicas made from phantom materials. Although the mechanics of brain tissue has been extensively studied, the selection of suitable phantom materials still relies on simplified metrics that fail to capture the highly nonlinear and rate-dependent behavior of the brain. In this work, we characterize the mechanical response at low-to-medium strain rates of common synthetic brain phantom materials, including Sylgard 527, Ecoflex Gel, and a composite hydrogel, and establish rigorous comparisons with porcine brain tissue as a benchmark. To this end, we employ a multi-step stress relaxation protocol under both uniaxial extension and compression across two different strain rates, introduce thermodynamic metrics, and perform comparisons based on rigorous statistical analysis. In particular, we isolate the thermodynamic equilibrium stress-strain curves, representing the baseline material's response, from rate-dependent dissipative inelastic effects, relatable to viscoelastic and poroelastic mechanisms. Our results provide a thorough comparison, revealing critical limitations in existing simulant materials and highlighting significant space for developing improved brain surrogates. The comprehensive dataset collected can also support the calibration of constitutive models, enhancing the predictive accuracy of numerical simulations involving brain tissue and its surrogates. More broadly, this study highlights the need to move beyond simplistic, scalar-metric-based comparisons of mechanical responses across materials.

Objective: To evaluate quantitative gait analysis in a sheep model (Ovis aries) of in-utero spina bifida repair. Method: Using a 3D marker-based motion capture system, we quantified spatiotemporal gait parameters and hindlimb joint angles in a sheep model of in-utero spina bifida repair and normal controls. A linear support vector machine classifier, trained and tested on subject-disjoint partitions with individual gait cycles as the unit of analysis, evaluated the relationship between hindlimb joint angles and clinician-assigned Texas Spinal Cord Injury Scale scores. Linear mixed-effects models compared spatiotemporal parameters, hoof height, and hip drop across 3- to 6-month evaluations and between groups. Results: In treated (n=14) and control (n=3) lambs, front- and hindlimb gait speed differed by timepoint (both p=0.025). Across 1,632 left and 1,661 right hindlimb gait cycles, held-out classification accuracy was 85.9% and 88.6% (Cohen's {kappa}=0.78 and 0.84) for left and right hindlimbs, respectively, against majority-class proportions of 76.5% and 75%. Left hindlimb hip drop differed between treated lambs and controls ({Delta}=0.8 cm, p=0.027). Between-group contrasts were considered exploratory given the small control group. Conclusion: Quantitative gait analysis was feasible in an in-utero spina bifida repair model. Hindlimb kinematic profiles classified clinical ambulation scores above the majority-class baseline.

The spring-loaded inverted pendulum (SLIP) model has effectively explained center of mass (CoM) and ground reaction force (GRF) dynamics during human gait despite its simplicity, describing both walking and running within the same mechanical principles. A recent study extended the explanatory capability of the SLIP model to joint dynamics by reproducing ankle joint dynamics during walking. However, during running, the model did not capture the empirical ankle torque profile, suggesting that such mechanical unification may not hold at the joint level. Because the previous model represented the body as a single point mass, it neglected leg inertia, which may become important during running as rapid motion generates substantial angular momentum. To address this limitation, we propose an extended two-mass model that separates the whole-body mass into the swing-leg and the remaining body mass, thereby explicitly incorporating swing-leg inertia into stance-leg joint dynamics. The swing leg is modeled using a linear spring and a torsional spring at the hip joint for its radial and rotational motion. Through inertial coupling, swing-leg dynamics influence the stance-leg ankle torque. Consequently, the proposed model reproduces the symmetric, single-peaked ankle torque profile observed experimentally while preserving the CoM and GRF of the previous model. These findings suggest that the discrepancy of the previous model during running may be explained by the omission of swing-leg inertia. By accounting for the gait-mode-dependent contribution of swing-leg inertia, the proposed model provides a unified mechanical framework for extending simple gait models from CoM-level behavior to joint dynamics.

Classical neuronal cable theory relies on quasi-static electric field approximations and neglects magnetic induction, Lorentz force coupling, and transient electromagnetic currents, limiting its ability to fully characterize action potential propagation within geometrically branched axons and dendrites. This work develops a coupled Maxwell-electromagnetic cable framework by integrating finite-difference time-domain (FDTD) solutions of Maxwells equations with extended Hodgkin-Huxley and Fitzhugh-Nagumo membrane dynamics, incorporating magnetic gating perturbations, electromagnetic trans-membrane currents IEM, and nanoscale quantum corrections for thin neural segments. Controlled propagation experiments are designed to quantify deviations from standard cable predictions across asymmetric and symmetric axonal bifurcation geometries. Numerical results demonstrate that inductive magnetic effects lower the critical branch radius for junction conduction failure and break symmetric action potential invasion in geometrically identical child branches under external transverse magnetic fields. An electromagnetic corrected geometric ratio GREM is proposed to revise impedance-matching conditions at branch points, accounting for size-dependent axial current imbalance induced by magnetic and displacement currents. Parent axon conduction velocity deviates substantially from the canonical [Formula] scaling law when electromagnetic feedback and quantum charge distributions are included, triggering early signal blockage at large cable diameters. Collectively, this study establishes that quasi-static cable models underestimate electromagnetic corrections to propagation speed, waveform shape, and bifurcation transmission fidelity; the coupled Maxwell-cable framework provides a comprehensive multi-physics tool for modeling electrodynamic signal behavior in complex neuronal architectures.

Bone strain and remodeling are important factors for bone health. Improved characterization of relationships between these factors in response to acute loading will improve our understanding of bone maintenance, adaptation to load, and pathophysiology. The purpose of this study was to evaluate colocalization between finite-element-estimated bone strain and [18F]sodium fluoride positron emission tomography imaging-estimated subchondral bone remodeling in the tibial diaphysis in response to acute exercise. We hypothesized that the exercise-induced change in standardized uptake values would be greatest in regions of higher strain. 26 participants underwent bilateral positron emission tomography/magnetic resonance scans before and after a bout of stair ascent and descent exercise. Finite element analysis was performed to calculate principal strain distributions in the tibial diaphysis. [18F]sodium fluoride standardized uptake values were calculated from pre- and post-exercise positron emission tomography data and mapped to the tibial surface. Mean pre-exercise and exercise-induced change in standardized uptake values were calculated and compared between regions defined by strain percentile. Both normalized mean pre-exercise and exercise-induced change in standardized uptake values increased with maximum principal strain in a pattern suggesting an exponential relationship. Pre-exercise and exercise-induced change in standardized uptake value were notably much greater in the region of highest maximum principal strains (estimated marginal means = 76.9 and 71.0 percentile, respectively) when compared to all other regions of the tibial diaphysis ([&le;] 61.1 and [&le;] 57.0 percentile, respectively, p < 0.001). This study is the first to directly show a relationship between bone mechanical loads and acute [18F]sodium fluoride uptake in humans in vivo. Our results highlight the influence of tensile strain magnitude on bone metabolic activity in response to acute loading. Furthermore, we demonstrate the utility of a combined finite element analysis and positron emission tomography approach to interrogate relationships between bone loading and bone metabolic response to acute exercise.

Computational models of the middle ear are typically developed within the linear range of the middle ear's response, which corresponds to typical sounds in daily life. Understanding middle-ear biomechanics under high-intensity sounds is also important for assessing potential auditory damage under adverse conditions. In this study, a finite-element (FE) model of the human middle ear plus the external ear canal is presented, employing a first-order Ogden hyperelastic material formulation and time-domain nonlinear analysis to capture nonlinear deformation of soft tissues and the resulting middle-ear response. The model was analyzed using pure-tone pressure excitations applied at the ear canal entrance over the 80 Hz to 8 kHz frequency range and 110 to 180 dB sound pressure levels (SPL). The ear canal pressure near the tympanic membrane (TM) and full-field three-dimensional response of the middle ear structures were examined. Across the investigated frequency range, the simulated middle-ear motion exhibited a transition from predominantly piston-like behavior at low frequencies, up to 1 kHz, to increasingly intricate multi-modal and multi-directional vibration patterns at higher frequencies. Nonlinear behaviors became increasingly significant at high frequencies and high stimulus levels, including compressive/expansive pressure-displacement relationships, waveform distortions, and the emergence of non-harmonic spectral components. The model developed in this study provides a computational framework for studying nonlinear sound transmission through the middle ear and its implications for high-intensity acoustic exposure. The findings contribute to a better understanding of middle-ear mechanics under adverse conditions and support future investigations of acoustic trauma and hearing protection strategies.